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ABSTRACT 

We have used the Green Bank Telescope (GBT) and Berkeley-Illinois-Maryland Association 
(BIMA) array to search for redshifted millimetre absorption in a sample of damped Lyman- 
alpha absorption systems (DLAs). This brings the number of published systems searched from 
18 to 30. In 17 cases we reach 3tr limits of r ^ 0.1, which is a significant improvement over 
the previous searches and more than sufficient to detect the 4 known redshifted millimetre 
absorbers (r > 1). While the CO rotational (millimetre) column density limits obtained are 
weaker than the electronic (optical) limits, they may provide useful limits below the atmo- 
spheric cut-off for the Lyman and Werner H2-bands in the UV (z a b s ~ 1.8). Using a model 
for the DLA metallicity evolution combined with assumed HCO + /H2 and CO/H2 conversion 
ratios, we use the molecular column density limits to calculate plausible H2 molecular frac- 
tion limits. Finally, we use these results to discuss the feasibility of detecting rotational CO 
transitions in DLAs with the next generation of large radio telescopes. 

Key words: quasars: absorption lines-galaxies: ISM-radio continuum: galaxies-cosmology: 
early universe 



1 INTRODUCTION 

Molecular absorption lines trace, and provide detailed physical and 
chemical information about, the cold dense component of the in- 
terstellar medium (ISM). Although many detailed studies exist for 
molecular clouds within our own Galaxy, only relatively recently 
has detailed information emerged for molecular abundances at 
high redshift, through absorption studies of redshifted UV molec- 
ular hydrogen lines (e.g. iLedoux et alJl2003tlReimers et al]|2003l) 
and mi llimetre-band rotational lines from molecular tracers (e.g. 
IWiklind & CombeJl996dY 

Observations of a range of different molecular transitions in 
gas clouds at high redshift would provide a wealth of informa- 
tion on star formation activity in external galaxies, potentially 
viewed at epochs when chemical abundances and environments 
were markedly different to today. Such information is invaluable 
for a detailed understanding of galactic formation and evolution. 
Furthermore, the narrowness of molecular lines reveals information 
about small-scale structure in the ISM. Comparison of relative line 
strengths yields information about the excitation mechanism, and in 
particular, probes the temperature of the Cosmic Microwave Back- 



ground (CMB), and hence the expected (1 + z)-dependence (e.g. 
Wiklind & Combes 1996c). Finally, comparisons of the relative ob- 
served frequencies of millimetre molecular lines (with each other, 
and/or with atomic transitions arising in the same cloud) can be 
used to check on any possible variation in certain combinat ions of 
the fu ndamental constants iDrinkwater et alj|l998t iMurphv et alJ 
2001). This last point was the prime motivation for our search for 
new redshifted millimetre absorbers reported in this paper, although 
we use the upper limits obtained to yield molecular fraction limits 
in low redshift DLAs, which is not possible using the UV H2 lines, 
which fall below the atmospheric cut-off of z a b s ~ 1.8. 

Curren tly, only 4 redshifted millimetre absorption systems are 
known (see Wiklind & Combeall999l and references therein), of 
which the highest redshift is 0.886 (PKS 1830-211). As a means 
of approaching a search for new high redshift radio absorbers 
system atically, we produced a catalogue of DLAs ICurran et alJ 
2002b) 1 , where large column densities (Nm > 2 x 10 20 cm" 2 ) 
are known to exist, and shortlisted those which are illuminated by 
radio-loud quasars (i.e. those with a measured radio flux density 
> 0.1 ly). This yielded 60 DLAs and sub-DLAs occulting radio- 



E-mail: sjc@phys.unsw.edu.au 



1 A version of this catalogue is kept updated on-line and is available from 
http://www.phys.unsw.edu.au/~sjc/dla 
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loud qu asars. Of these, 37 have be en searched for 21 -cm absorp- 
tion (see Kaneka rChengaludl2003HCurran et all2 004b). Selecting 
those of 12-mm and 3-mm flux densities > 0.1 Jy, gives 18 sys- 
tems wMch_have previously been searched for millimetre absorp- 
tion ICurran et all2003l) . this number now being increased in total 
to 30 with the observations we present in this paper. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 12-mm observations 

In Tabled we list all of the known DLAs and sub-DLAs which are 
illuminated by background sources of ^ 0.1 Jy at centimetre wave- 
lengths and have molecular transitions which could be detected in 
the GBT K-band. In all cases this is both the HCN and HCO+ 
J = — > 1 rotational transitions' 5 . We chose to observe the latter 
transition since this tends to be slightly stronger than HCN in the 4 
known millimetre absorption systems. 

Apart from the sources listed in TableQ there are only 3 other 
such DLAs known which have millimetre transitions redshifted to 
« 12-mm: Q 0438-436, RX J1028.6-0844 andLBQS 1213+0922. 
Although strong (S ~ 3 Jy at 12-mm), the first source is far too 
south to observe with the GBT, and the latter two were estimated 
to have flux densities too low to yield good optical depth limits 
within a reasonable integration time. Therefore, using the GBT we 
hoped for a near complete survey of 12-mm absorption in known 
DLAs, although due to time constraints we did not observe some of 
the sub-DLAs <Turnshek et alll989l : lLanzetta et alll99ll) towards 
1017+109, 1021-006 and 1402+044, nor the N m = 4 x 10 20 cm " 2 
absorber at z abs = 3.533 towards 1251-407 <EUison et all200lh . 

The observations were performed in April 2003, during sev- 
eral consecutive days. Typically, under acceptable observing condi- 
tions, the system temperature was 40-100 K. At the time of observ- 
ing the wide bandwidth spectrometer was unfortunately unavailable 
and so we used the digital spectral processor with an observing 
bandwidth of 40 MHz, over 1024 channels. This gave a redshift 
coverage of Az a b s ~ ±0.004, sufficiently wide to cover the ab- 
sorption redshift uncertainty in most of the DLAs. However, the 
optical redshifts of 1614+051 and2131-045 are poorly constrained 
(Table H) and so we observed these with five and three overlap- 
ping bands, respectively. Dual polarization was used to optimise 
the sensitivity and beam switching facilitated background emission 
removal. 

The data were reduced using the AIPS++ single dish package 
and baselines were removed from the spectra by fitting a polyno- 
mial function. The continuum levels were measured either during 
GBT pointing scans of the brightest sources or with the Very Large 
Array (VLA) in its most compact (D) configuration. These latter 
data were taken in May 2003, and thus the continuum levels might 
have altered due to quasar variability, although, as seen from Ta- 
ble [3] these appear to be stable over the 9-11 month period since 
the ATCA observations. For sources in which no continuum levels 
were successfully measured within a month of the GBT observa- 
tions, flux densities were taken or estimated from the literature (see 
Table|3j. 

2 Unless otherwise stated, B 1950.0 names are used throughout this paper. 

3 Due to the high dipole moment of the density tracing molecules (i.e. HCN 
and HCO+), these will generally be close to the ground state causing the 
opacity to rival that of the traditional molecular tracer CO which will be 
excited to higher rotational levels by the CMB at these redshifts. 



Table 1. All of the 22 GHz illuminated DLAs and sub-DLAs for 
which a strong molecular rotational transition falls into the G BT K- 
band (range 18-26.5 GHz): Prior to the results of iRantakvro et all 
|2004|) . all of the flux densities were measured or estimated to be 
j5 0.1 Jy at 12-mm. Nm is the neutral atomic hydrogen column den- 
sity (cm -2 ) and z a b s is the redshift of the DLA, respectively. u a ^ s is 
the observed frequency of the HCO+ — » 1 line and S is the flux 
density of the background quasar previously measured at 12-mm. 
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(a) AT CA calibrator, (6 ) Kovalev et alJ <1999l). (el iTerasranta et all 
Il998h . Idl lCurran etall <2003h . (el lRantakvro et al]|2004l) . Note that 
som e flu x densities* were estimated from neighbouring frequencies 
(see ICurran et alj|2002tl) and that 1614+051 is a Gigahertz Peaked 
Spectrum source. 



Also in May 2003 HCN -f 1 at z abs = 2.713 towards 
1402+044 was observed with the H22 receiver on the Nobeyama 
45-m during less than ideal weather conditions during observations 
for a different (3-mm) project jCurran et all2004cj) . The 250 MHz 
AOS-W spectrometer (0.25 MHz channels) gave a redshift cover- 
age of Az abs ~ ±0.02. 

2.2 3-mm observations 

In addition to the GBT survey, follow-up CO observations of can- 
didate molecular absorption detections ICurran et alJl20023l were 
carried out using the BIMA array (Table|2j- The spectrometer was 
configured in mode 5, with a central high resolution 100 MHz 
wide window (0.78 MHz channels) flanked by two lower res- 
olution 200 MHz wide windows (3.125 MHz channels), giving 
Az abs « ±0.003 - 0.007 for z ~ 0.2 - 2. Some overlap was 
added in order to eliminate edge channels and ensure the proper 
matching of the windows. The antennas where tuned to the red- 
shifted frequency of the chosen CO transition (Lovas 1992), using 
a redshift that was the average of that of the reported candidate 
absorption and the corresponding HI absorption (Az a b s ^ 0.001 
in all cases). To correct for possible spurious absorption detections 
due to RF or IF passband features, we interleaved observations of a 
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Table 2. The BIMA array targets. The parameters are as 
given in Table [Tl where u ^ s is the observed frequency of 
the CO -> 1 line (2^3 for 0458-020) and S is the 
flux density of the background quasar previously measured 
at 3-mm. 
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passband calibrator (a strong unresolved radio source) with those of 
the target quasars every hour. The total integration time on the cali- 
brator is such that the noise introduced by the passband calibration 
is negligible. Since the data could be self-calibrated, no observa- 
tions of phase calibrators were performed. The data were reduced 
using the MIRIAD interferometry reduction package. The visibil- 
ity cubes were imaged using natural weighting corrected by system 
temperature, the resulting images were cleaned, and spectra of the 
central positions were obtained. The final passband-corrected spec- 
tra were produced by dividing the spectrum of the source by that of 
its calibrator, then restoring the continuum level by multiplying by 
the mean value of the calibrator. The resulting spectra, along with 
the GBT results, can be viewed at synergy url here please 
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Figure 1. The normalised CO column density versus the absorption red- 
shift. The unfilled markers designate upper limits - squares for the mil- 
limetre searches and circles for the optical searches. The filled stars repre- 
sent the 4 known millimetre absorption systems. The lookback time is for 
fl m = 0.30, tt A = 0.70, H = TOkms" 1 Mpc" 1 . 



line- width (FWHM) of 10 km s -1 , which is close to those of the 
4 known systems (e.g. IWiklind & CombesllT996cl) . This is done 
for a resolution of the same value and so the column density 
limits, which are calculated for an excitation temperature 7 of 10 
K at z a bs = 0, represent a one channel 3a detection of a 10 
km s _1 wide line. 



3 RESULTS AND DISCUSSION 
3.1 Search results 

In Table|3|we show the best previously published optical depths to- 
gether with our new results, for which, despite the improved limits, 
there are no absorption features of ^ 3<r over ^ 3 km s~\ the full 
resolution of the BIMA array observations. For all of the limits we 
estimate the 3<r upper limits on the total column density of each 
molecule from 

where v is the rest frequency of the J — > J + 1 transition, gj+i and 
Aj+i-tj are the statistical weight and the Einstein A-coefficient 
of the transition, respectively, Q = Y2'j'=o9 J e~ Ej ^ kTx is the 
partition function 5 , for the excitation temperature, T x , and J rdv 
is the 3cr upper limit of the velocity integrated optical depth of the 
line 6 . 

Since the derived optical depth limits depend up on the r.m .s. 
noise and thus the velocity resolution (see Curran et al. 20 02al) . in 
our previous articles we presented all of the optical depth limits 
normalised to a spectral resolution of 1 km s~\ the finest resolu- 
tion typical of most current wide-band spectrometers. Previously, 
we also quoted column density upper limits per unit km s _1 line- 
width. However, in order to use our limits to constrain molecu- 
lar fractions at low redshift (Section 3.2), we shall now adopt a 

4 These are taken from Chan dra et all jl995lll996l) or derived from the 
dipole moment (e.g. lRohlfs & Wilsorl200d) . 

5 The e nergy of each leve l, Ej, is obtained from the JPL Spectral Line 
Catalog iPickett et alll998l) . 

6 Jrdv « 1.06 7"mm x FWHM for a Gaussian profile, where T mm is 
the 3<r peak optical depth limit (Tablel3*l. 



3.2 Deriving molecular fraction limits 

Despite our improved limits there has yet to be a detection of a ro- 
tational molecular transition in a DLA. In many cases the limits far 
exceed those required to detect the 4 known redshifted millimetre 
absorbers (see below, Figs.Qand|2j. These systems have molecular 
fractions of / = 2jVH + 2 Nhi ~ 0.3 — 1.0 (e.g. lCombes & Wiklind 

1998), c.f. / ~ 10~ 7 - 10" 2 for the 10 DLAs in which molec- 
ular hydro gen has b een detected at UV and optical wavelengths 
( Reimers et al. 2003 and references therein). In this section we in- 
vestigate the possibility of detecting molecular tracers in systems 
with such low molecular fraction. 

In F ig. [T] we show t he CO rotational (Table 13 and electronic 
iGe et alJll997t iLu et alj|l999l IPetitiean et alj|2002l) column den- 
sity limits along with the values for the 4 known systems. While this 
illustrates the sensitivity to detecting absorption in the known sys- 
tems, it clearly shows that the radio searches are considerably less 
sensitive than the optical CO searches 8 . They do, however, comple- 
ment the optical results in giving limits at redshifts of z a b s ~ 1.8 
(below the atmospheric cut-off for direct optical H2 detections). In 
Fig.|2|we show the HCO + column density limits. Again we see that 
our limit s are more than sufficient to detect HCO + in the 4 known 
systems JWiklind & CombeJl99lll996albiri997l) . with two cases 
having been searched to at least an order of magnitude better than 
that required to detect the known systems at low redshift. 

7 Due to the increase of CMB temperature with redshift, T^mb = 
2.73(1 + 2; a bs)> Tx increases to 20 K at z a j, s = 3.75, the highest red- 
shift of our sample. The main effect of this is to increase the HCO+ — > 1 
column density estimates by a factor of ^ 4 in comparison to a constant 
value of 10 K. 

8 Although the deep integrations of 0235+164 and 0458-020 with the NRO 
45-m telescope iTakahara et al. 1987; Tsuboi & Nakai 1991) do approach 
the optical limits. 
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Table 3. Summary of published searches for molecular tracer absorption in DLAs and sub-DLAs. i/ obs is the approximate observed frequency 
(GHz), V is the visual magnitude of the background quasar and S is the flux density (Jy) at u ^ s which, unless flagged, is the value obtained 
during the actual observations (the limits are Icr and blanks in this field indicate that no value is given in the literature). Ahi (cm~ 2 ) is the DLA 
column density from the Lyman-alpha line and T21 cm is the normalised peak optical depth of the redshifted 21-cm line (see lCurran et all2 004b 
for details). The optical depth of the relevant millimetre line is calculated from r = — ln(l — 3<r rms /5), where cr TInB is the r.m.s. noise level. 
Since cr rms is dependent on the spectral resolution, we take the various published values and recalculate (r rms at a resolution of 10 km s~ 1 (r mm ) 
and quote only the best existing limit. For all optical depths, 3cr upper limits are quoted and "— " designates where 3<r > Scout, thus not giving 
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19.2 


« 0.07 A ' V 


3 x 10 20 


< 0.015 


< 0.08 


< 1 x 10 13 


10 




2.996 






22.5 










< 0.1 


< 2 x 10 13 


10 


1402+044 


2.713 


HCO+ 


-> 1 


24.0 


19.8 


0.21* 


8 x 10 19 




< 0.008 


< 1 x 10 12 


10 






HCN0 - 


-> 1 


23.9 




0.6 






< 0.02 


< 8 x 10 12 


10 


1418-064 


3.449 


HCO+ 


-> 1 


20.1 


18.5 


0.14 


3 x 10 20 




< 0.02 


< 4 x 10 12 


10 


1451-375 


0.2761 


HCO+ 1 


-> 2 


139.8 


16.7 


0.6 


1 x 10 20 


< 0.007 


< 0.2 


< 7 x 10 12 


8 






COO - 


1 


90.3 




1.2 






< 0.1 


< 8 x 10 15 


8 
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Flux densities: A ATCA June & August 2002 iCurran et all 2003: Rant akvro et all2004l) . V VLA May 2003. Where the flux densities could not 
be obtained from the data/article - a the BIMA array value, b interpolated between 1 1 GHz and K-band, interpolated between radio and 0.5-10 
keV, d extrapolated from 0.4 and 5 GHz. "This observation of 1402+044 is from preliminary (February 2003) GBT observations using the wide 
band spectrometer. 
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Table 3. Continued 



DLA 


2 b 


Transition 


Vobs 


V 


S 


AT H I 


T21 cm 




Nmm 


Ref. 


1614+051 


2.52 


HCO+ 


-> 1 


25.3 


19.5 


~ 0.3 C 


3 x 10 20 




< 0.03 


< 4 x 10 12 


10 


2131-045 


3.27 


HCO+ 


-» 1 


20.9 


20.0 


< 0.06- 4 


1 x 10 20 








10 


2128-123 


0.4298 


CS 2 - 


> 3 


102.8 


15.5 


0.75 


2 x 10 19 


< 0.003 


< 0.04 


< 1 x 10 13 


7 






CS 3 - 


> 4 


137.1 




0.70 






< 0.05 


< 2 x 10 13 


7 






CO 1 - 


♦ 2 


161.2 




0.5 






< 0.2 


< 6 x 10 15 


7 






CO 2 - 


♦ 3 


241.9 




0.4 






< 0.2 


< 1 x 10 16 


7 


2136+141 


2.1346 


HCO+ 2 


-* 3 


85.4 


18.9 


0.59 


6 x 10 19 




< 0.1 


< 4 x 10 12 


6 






CO 2 - 


♦ 3 


110.3 




0.27 






< 0.2 


< 9 x 10 15 


7 






CO 3 - 


* 4 


147.1 




0.21 






< 0.2 


< 2 x 10 16 


7 






CO 5 - 


♦ 6 


220.6 




0.16 






< 0.6 


< 8 x 10 17 


7 


2342+342 


2.908 


HCO+ 


-> 1 


22.8 


18.4 


~ 0.1 d 


2 x 10 21 


< 0.03 


< 0.04 


< 7 x 10 12 


to 
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Figure 2. The normalised HCO + column density versus the absorption red- 
shift. Again the stars represent the 4 known systems and the inverted trian- 
gles are the DLA upper limits (Tablel3l. 

Although we cannot match the sensitivity to molecular ab- 
sorption provided by the optical results, we can nevertheless at- 
tempt to derive constraints on the molecular fraction for absorption 
systems at low redshift. The major obstacle in this is determining 
which A r co-A r H 2 conversion ratio to apply: Unlike HCO + (dis- 
cussed l ater), the ratio may va ry from A/h 2 ~ 10 6 A r co for dif- 
fuse gas iLiszt & Luca s 2000a) to Nn 2 ~ 10 4 iVco for dense, dark 
Galactic clouds. This latter value is applied to the 4 known millime- 
tre absorbers to yield molecular fractions from rotational lines (see 
Wikli nd & Combe JT 998 ) 9 . However, DLAs have lower metallici- 
ties than Galactic systems and are less visually obscured than the 4 
known systems, thus casting doubt on whether applying this con- 
version ratio is justified. 

If we extrapolate the metallicity ([M/H] 10 ) evolution of the 
general DLA population 1 according to [M /Hl ~ — 0.26z a b s 
jKulkarni & Fall2002t|Prochaska et all2003l) . we obtain [M/H] « 

9 Since all of these sources absorb at low redshift, z a t, s $C 0.87, and have 
V £ 20 there is currently little chance of detecting the H2 line directly. 

10 Defined as the heavy element abundance with respect to hydrogen, rela- 
tive to that of the solar neighbourhood: [M/H] = log 10 [iV(M)/JV(H)] - 
log 10 [iV(M)/7V(H)] Q . " 

11 Due to many of the sources in Tablel3lhaving low redshifts and/or high 
visual magnitudes, there are only metallicity measurem ents for half of these 
sources, all of which belong to the sample of Prochaska et al. 1 2003). There- 
fore, by definition these follow the same metallicity evolution of the gen- 
eral DLA population rather than that of the H2-bearing DLAs 1 Curran et al. 
l2004dl . 



Lookback time [Gyr] 
6 7 8 9 10 11 11.5 




I 1 1 I 1 1 1 1 I 1 1 1 1 1 1 i_l 

1 2 3 

Absorption redshift 

Figure 3. The metallicity evolution of DLAs. The filled circles represent 
the DLAs which ex hibit H2 absorption with the least-squares fit shown 
ICurran et af]|2004d) . the small unfilled diamonds represent the metallic- 
ity measurements of 60 DLAs betwee n 1.9 < z a .hn < 3.4 (i. e. the general 
population), again with the fit shown 1 Prochaska et al. 2003), and the un- 
filled squares are the DLAs searched for millimetre absorption (Tablel3*l. 



-0.75 at ~ (Fig.0. At zero redshift this is still significantly 
lower than solar values suggesting that, in addition to lower enrich- 
ment at high redshift, there may be appreciable dust depletion of 
the m etals, as in dense Galactic clouds (e.g. ISavage & Se mbach 
1996), although this may be more rele vant to the H2 -bearin g DLAs 
wher e / > 10~ 2 at_Za ba < 1.8 JCurran et al]l2004ch . More- 
over, iLedoux et alJ 120031) suggest that the large depletion factors 
in DLAs indicate a significant dust content. Therefore, a fair com- 
promise between the diffuse/dark cases may be to assume the con- 
version ratio for dark clouds but scaled by the abundance of metals 
from which the tracer molecules form, i.e. 7Vh 2 ~ 10 4 ' 75 A r co at 
-Zabs ~ 0. Since the metallicity decreases with increasing redshift, 
we may expect the conversion of tracer to molecular hydrogen col- 
umn density to scale accordingly. We therefore apply the ratio of 
A*h 2 ~ 10 0,26ZabB+4,75 A i 'co, based upon the metallicity evolution 
of the general DLA population 12 . This gives 7Vh 2 ~ 10 5 A r co at 
z ~ 1 rising to A^H q ~ 10 6 A r co at z ~ 5. F rom optical CO limits, 
iBlack et alj<1987l) and lChaffee et alJll988h have previously noted 



Since the molecular hydrogen fraction shows a strong ant i-correlation 
with redshift for the H 2 -bearing DLAs iCurran et aT] l2004c and Fig. 
we may expect less depletion and thus relatively higher metallicities. Fig. [3] 
does appear to suggest, however, that the metallicity is dominated by poorer 
chemical enrichment at higher redshifts. 
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Figure 4. Molecular fractions measured in redshifted absorbers. Again the 
unfilled symbols represent the upper limits from the millimetre searches, 
Nh 2 = 10 4 75 A r co and T x = 10 K at z abs = 0, with the squares 
representing CO and the inverted triangles HCO+. The circles represent the 
H2 -bearing DLAs and the stars the 4 known millimetre absorption systems 
(Nn 2 = 10 4 N CO and T x = 10 K, Vz abs ). A log version of the plot is 
also shown by means of a detail of the H2 -bearing systems. This clearly 
illustrates the evolution of / described bv lCurran etaT]|2004j| . 



that the carbon to hydrogen column density ratio is a tenth of the 
Galactic value at z abs = 1.7 and 2.3. 

In addition to CO, we can use the HCO + limits, which 
could be a better choice of tracer, since a near constant ratio of 
jV HCO+ = 2 - 3 x 10~^Ho is found over various regimes in 
the Galaxy iLiszt&Lucasl EoOObl. Applying the conversion thus 



suggests that 7Vh 2 



10 M 



lO^iV, 



HCO- 



at z 



5, i.e. 



7Vh 2 » 10 ' 26z '» b3+9 - 34 A r HC o+. In Fig. @ we show the derived 
molecular fractions of the CO and HCO + limits together with those 
of the H2-bearing DLAs 13 . From this we see that the high con- 
version ratio and excitation above the — ► 1 transition make the 
HCO + molecule insensitive with current telescopes to / < 1: Al- 
though in many cases our GBT observations yield significantly bet- 
ter limits than previously (Table[3j, these are all at high redshift and 
therefore convert less favourably to molecular fractions due to the 
steep metallicity evolution. Thus, we obtain near constant limits to 
the fraction across the entire redshift range searched for HCO + . 
Note that the optical CO limits (Fig.0 occupy the space between 



the HCO + limits and the H 2 -bearing DLAs (1.7 



4.4, 



—3 ~ log 10 / 5; —0.5) when the conversion is applied. 

We do see however, based on our redshift dependent conver- 
sion factor, that the millimetre CO searches generally give interest- 
ing (low redshift) limits, with two of the limits within the range 
of the H2 -bearing DLAs. This suggests that, were these DLAs 
H2-bearing 14 , making the Nqo-Nyi 2 conversion reasonable, the 



13 Versions of Fig.|4]for A r jj 2 = 10 6 Af C o and 10 4 AT C O can be viewed 
at synergy url here please 

14 H2 is detected in S; 20% of recent DLA surveys <Ledoux et all2003h . 



CO rotational transition could have been detected in these DLAs: 
0235+164 (z abs = 0.524) and 0738+313 (2 abs = 0.221). Un- 
fortunately, ther e are no metallicity m easurements available for 
these, although (Prochaska et al. 2003) use the value of [M/H] = 
—0.22 ± 0.15 for 0235+164 whi ch comes from model dependent 
X-ray measurements fTurnshek etafl[ 2003 ) . Several other models 
give lower values than this, thus suggesting that this DLA may 
belong to the general population ( [M/H ] « —0.9), although the 
value quoted by Prochaska et al. (2003) could place it in either 
population (Fig. The general population fit is, however, con- 
sistent with thejan^e_rrf-1.4_<_rM/H] < -0.4 at z abs ~ 0.4 
estimated by Tumsh ek et alj 120031) and the fact that we did not 
detect CO. Regarding the molecular fraction, the data are sensi- 
tive to log 10 / ~ —1.5 for the general population metallicity and 
log 10 / ~ —2.1 for the value used bv lProchaska et alj (;2003). This 
latter fraction shifts this DLA further into the Lb-bearing DLA 
regime. 

For more diffuse clouds of low dust content our limits become 
weaker with / ^ 0.4, c.f. ^ 0.03 previously, being obtained for 
JVh 2 ~ 10 6 iVco and T x = 10 K at z abs = 0, although in the 
diffuse gas case the column density limits would improve slightly 
due to lower excitation temperatures (at best a factor of ~ 0.5 lower 
than given in Table|3j 13 . That is, the clouds being diffuse could also 
explain why no millimetre absorption of CO was detected. 



3.3 Detecting millimetre lines in DLAs with the next 
generation of large radio telescopes 

Applying the conversion of tracer column densities described 
above, we can ascertain the likelihood of detecting these molecular 
tracers in DLAs with the next generation of large radio telescopes. 
For iV H2 = 10 4 ' 75 iVco and T x = 10 K at z abs = 0: 

(i) The Atacama Large Millimeter Array (ALMA) 15 : The 3-mm 
band of this telescope will give CO — > 1 coverage to z ^ 0.34. 
At 10 km s _1 resolution, one hour of integration will give an r.m.s. 
of ~ 0.3 mJy. For a 0.3 Jy continuum flux (fairly typical in Table 
[5), this corresponds to a 3<r column density limit of iVco ~ 10 14 
cm -2 per channel. Applying the conversion at these redshifts gives 
JVh 2 ~ 10 19 cm" 2 thus giving / ~ 10" 3 - 0.1 for DLAs. If we 
extrapolate the molecular fractions in H2 -bearing systems back to 
z 1z 0.3 (Fig.[4j, we see that such limits could well be sufficient to 
detect rotational CO in DLAs at low redshift. 

(ii) The Square Kilometre Array (SKA) 16 : With a tuning range 
of 0.15 to 20 GHz this telescope can observe CO at z ^ 4.8. With 
the same parameters as above, an r.m.s. of 5 /iJy is reached giving 
Afco ~ 10 12 cm -2 per channel 17 , thus rivalling current optical 
limits. At z ~ 5, 7Vh 2 ~ 10 6 iVco giving / ~ 10~ 4 - 10~ 2 for 
DLAs, i.e. the molecular fraction range of 7 of the 10 H 2 -bearing 
DLAs (e.g. iReimers et al] |2003). Unfortunately at z Z 5, extrap- 
olating the molecular fraction-redshift anti-correlation, we expect 
/ ^ 10 -12 , thus making it unlikely that redshifted millimetre tran- 



and so we perhaps expect a rj 1/5 chance of having observed an H2- 
bearing system. 

1 5 http ://w ww. eso.org/proj ects/alma/ 

16 http://www.skatelescope.org/ 

17 Again, this is for a flux density of 0.3 Jy. Due to the antenna temperature 
dominating the system at these low frequencies, an increased flux does not 
improve the optical depth to the same extent as for the millimetre observa- 
tions. Therefore, the column density estimate is fairly robust. 
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sitions could be detected with the SKA, despite its superior sen- 
sitivity 19 . 



4 SUMMARY 

We have used the GBT and BIMA array to search for HCO + and 
CO absorption in 19 damped Lyman-alpha absorption systems. Al- 
though in many cases the sensitivity reached far surpasses that re- 
quired to detect the 4 known redshifted millimetre absorption sys- 
tems, a rotational millimetre-band transition has yet to be detected 
in a DLA. From optical studies it is revealed that molecular frac- 
tions in DLAs are much lower than those of the 4 known systems 
and, in those where H2 has been detected, the H2 molecular frac- 
tion in DLAs decreases with increasing redshift. 

This evolution is also apparent in DLA metallicities, which 
we use to gauge the abundance of heavy elements available to form 
tracer molecules. By scaling the H2-CO & -HCO + column density 
ratios accordingly, we find: 

• Although the GBT results offer a significant improvement 
over previous HCO + searches, the high redshifts (z > 2.5) yield 
correspondingly large conversion ratios and excitation tempera- 
tures, i.e. Nh 2 ~ 10 10 iV HCO + , rendering high redshift searches 
sensitive only to molecular fractions of / ~ 1. 

• Since many of the CO searches are performed in the 3-mm 
band, these are generally at much lower redshift. So in addition 
to the z jj 1.7 optical limits of electronic transitions, we can 
estimate lower redshift values from the rotational transitions. Al- 
though the statistics are small and the quality of the limits var- 
ied (i.e. some limits are only sensitive to / Z 0.9), assuming 
JVh 2 w 10°' 26ZabB+4 ' 75 iVco, we obtain, for the general (non-H 2 ) 
DLA population, a 3u limit of / ~ 0.03 at z = 0.5. For more dif- 
fuse clouds (JVh 2 ~ 10°- 26Zaba+6 iVco) this increases to / < 0.2. 

Finally, we discuss the possibility of detecting CO absorption with 
the next generation of large radio telescopes in the context of the 
metallicity evolution model. Although the SKA will reach sensitiv- 
ities of ~ 5 /iJy after one hour of integration (c.f. ~ 300 /iJy with 
ALMA), ALMA is the most likely to detect millimetre absorption 
in DLAs at low redshift due to its high frequency capabilities. 
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